91 research outputs found
Numerical solution of the unsteady Navier-Stokes equation
The construction and the analysis of nonoscillatory shock capturing methods for the approximation of hyperbolic conservation laws are discussed. These schemes share many desirable properties with total variation diminishing schemes, but TVD schemes have at most first-order accuracy, in the sense of truncation error, at extrema of the solution. In this paper a uniformly second-order approximation is constructed, which is nonoscillatory in the sense that the number of extrema of the discrete solution is not increasing in time. This is achieved via a nonoscillatory piecewise linear reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell
Optimal Data Collection For Informative Rankings Expose Well-Connected Graphs
Given a graph where vertices represent alternatives and arcs represent
pairwise comparison data, the statistical ranking problem is to find a
potential function, defined on the vertices, such that the gradient of the
potential function agrees with the pairwise comparisons. Our goal in this paper
is to develop a method for collecting data for which the least squares
estimator for the ranking problem has maximal Fisher information. Our approach,
based on experimental design, is to view data collection as a bi-level
optimization problem where the inner problem is the ranking problem and the
outer problem is to identify data which maximizes the informativeness of the
ranking. Under certain assumptions, the data collection problem decouples,
reducing to a problem of finding multigraphs with large algebraic connectivity.
This reduction of the data collection problem to graph-theoretic questions is
one of the primary contributions of this work. As an application, we study the
Yahoo! Movie user rating dataset and demonstrate that the addition of a small
number of well-chosen pairwise comparisons can significantly increase the
Fisher informativeness of the ranking. As another application, we study the
2011-12 NCAA football schedule and propose schedules with the same number of
games which are significantly more informative. Using spectral clustering
methods to identify highly-connected communities within the division, we argue
that the NCAA could improve its notoriously poor rankings by simply scheduling
more out-of-conference games.Comment: 31 pages, 10 figures, 3 table
PDEs with Compressed Solutions
Sparsity plays a central role in recent developments in signal processing,
linear algebra, statistics, optimization, and other fields. In these
developments, sparsity is promoted through the addition of an norm (or
related quantity) as a constraint or penalty in a variational principle. We
apply this approach to partial differential equations that come from a
variational quantity, either by minimization (to obtain an elliptic PDE) or by
gradient flow (to obtain a parabolic PDE). Also, we show that some PDEs can be
rewritten in an form, such as the divisible sandpile problem and
signum-Gordon. Addition of an term in the variational principle leads to
a modified PDE where a subgradient term appears. It is known that modified PDEs
of this form will often have solutions with compact support, which corresponds
to the discrete solution being sparse. We show that this is advantageous
numerically through the use of efficient algorithms for solving based
problems.Comment: 21 pages, 15 figure
An L1 Penalty Method for General Obstacle Problems
We construct an efficient numerical scheme for solving obstacle problems in
divergence form. The numerical method is based on a reformulation of the
obstacle in terms of an L1-like penalty on the variational problem. The
reformulation is an exact regularizer in the sense that for large (but finite)
penalty parameter, we recover the exact solution. Our formulation is applied to
classical elliptic obstacle problems as well as some related free boundary
problems, for example the two-phase membrane problem and the Hele-Shaw model.
One advantage of the proposed method is that the free boundary inherent in the
obstacle problem arises naturally in our energy minimization without any need
for problem specific or complicated discretization. In addition, our scheme
also works for nonlinear variational inequalities arising from convex
minimization problems.Comment: 20 pages, 18 figure
Triangle based TVD schemes for hyperbolic conservation laws
A triangle based total variation diminishing (TVD) scheme for the numerical approximation of hyperbolic conservation laws in two space dimensions is constructed. The novelty of the scheme lies in the nature of the preprocessing of the cell averaged data, which is accomplished via a nearest neighbor linear interpolation followed by a slope limiting procedures. Two such limiting procedures are suggested. The resulting method is considerably more simple than other triangle based non-oscillatory approximations which, like this scheme, approximate the flux up to second order accuracy. Numerical results for linear advection and Burgers' equation are presented
- …